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ABSTRACT 


Experimental data for gas hydrate formation in the methane- 
isobutane-water system have been obtained over a range of gas phase 
compositions from 0.4 to 63.6 mole percent isobutane. The highest 
composition for which data previously existed was 4.6 percent isobutane. 
The range of experimental pressures was from 23.0 to 1,460 psia and 
the range of temperature was from 33.2 to 68.8°F. Experimental data 
were also obtained on the initial hydrate formation conditions for the 
water rich liquid-hydrocarbon rich liquid-hydrate equilibrium for 
three mixtures of light: hydrocarbons containing methane, ethane, 
propane and isobutane. 

A recently proposed prediction method based on the Kihara cell 
potential model has been used to predict gas hydrate formation conditions 
with mixtures of methane and isobutane. The predicted pressures agreed 
well with experimental values at high isobutane concentrations but were 
considerably higher than the experimental values at lower compositions. 
For the liquefied hydrocarbon mixtures, the predicted temperatures 
were within 2.8°F of the experimental values at the bubble point 


pressures. 
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CHAPTER 1 


INTRODUCTION 

Gas hydrates, resembling snow or loose ice in appearance, are 
solid crystalline substances having the general formula M.nH,0, where 
one or more hydrate forming molecules M associate with n "host" water 
molecules. The metastable water crystal lattice is Sabi ced by 
physical forces between the "host" and "guest" molecules. 

The important role of gas hydrates in gas processing industries 
was first realized in 1934 when Hammerschmidt [1] discovered that the 
hydrates of methane, ethane, propane and isobutane were responsible 
for the difficulties in the natural gas industry in the measurement 
and transportation of the gas. The movement of the gas through the 
pipeline tends to collect and compress the hydrates at low spots until 
the line may become entirely plugged. Research for the prevention of 
hydrate formation and, consequently, dehydration of the gas has been 
carried out since that time. Considerable effort has been expended 
in finding the conditions of temperature and pressure at which hydrate 
will form, as well as in determining the equilibrium moisture content 
data for the gas. There is additional interest in gas hydrate forma- 
tion connected with a series of processes which are carried out in 
other industries, for example, in distillation columns due to accumulat- 
ed moisture and in heavy water production processes involving mixtures 
of hydrogen sulfide and water. Several prospective applications for 
the industrial use of gas hydrate include (1) the desalination of sea 
water; (2) the storage of gases; (3) the separation of gas and liquid 


mixtures. Therefore, the prediction of initial hydrate formation 
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condition becomes important in relation to several industrial processes. 

There are many experimental data reported in the literature for 
gas hydrate formation conditions. By using these reported data several 
prediction methods have been developed and evaluated for their useful- 
ness. Although much data exists, that available for mixtures of 
methane and isobutane is very scanty. Therefore, one objective of 
this work was to extend and complete the hydrate formation data for 
methane-isobutane mixtures. 

Recent work, notably that of Parrish and Prausnitz [2], in 
developing a prediction method based on the statistical-mechanical] 
theory of hydrate formation has opened the way to more effectively 
predicting hydrate formation in gas-liquid systems and has made it 
possible to attempt the prediction of hydrate formation in liquid- 
liquid systems. Consequently, another objective of this work was to 
investigate the phase behavior of hydrate formation conditions for 
liquefied light hydrocarbons and to determine at what condition and 
for what kinds of components in the mixtures hydrate formation condi- 


tions can properly be predicted with the new method. 
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CHAPTER 2 


LITERATURE REVIEW 

Von Stackelberg and Muller [3,4], Claussen [5,6], Pauling and 
Marsh [7], and Jeffrey and McMullan [8] studied the crystal structures 
of gas hydrates and concluded that for most nonpolar and some weakly 
polar molecules gas hydrates crystallise in either of two structures, 
referred to as Structure I and Structure II. The physical properties 
of these structures are given in Table 1. The structure formed depends 
primarily on the size of the guest molecule. When the effective 
diameter of the hydrate-forming molecules is less than 5.1 A, all the 
cavities in hydrate Structure I can be filled and the ideal formula 


would be M,.5 3/4 H,0 (M, = Ar, Kr, CH S; etc.). If the effective 


ny 
° 
diameter of the hydrate-forming molecule is less than 5.8 A but greater 


than 5.1 A then the hydrate-forming molecules can fill only six large 


cavities and would have the ideal formula M, +7 CBS H,,0 (M, = SO 


a ae 
CoH» CoH, Clo, Br,» etc.). When two types of molecules are present 


Simultaneously in the gas mixture, a so called "mixed hydrate" with 


the ideal formula 2M, . 6M, .46H.0 will form. Similarly, if the effective 


diameter of the hydrate-forming molecules is greater than 5.8 A but 


4 Fe CHa » 1-Cy 10? 


II will form with the ideal formula M,-17H,0. For a mixture Mo + Mas 


less than 6.7 A (M H CHABr. » etc.) a hydrate Structure 
if the effective diameter of M. is less than 5.0 A, the "mixed 
hydrate" will form with the ideal formula 2M4.M,-17H,0. Both types of 


mixed hydrates are considered as solid solutions [9] in which the smaller 


molecules help to stabilize the hydrate structure. Any gas having a 
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molecular dimension larger than 6.7 A such as n-butane and larger 
molecules are too big to form hydrates. It has been found that the 
quantum-gas molecules (helium, hydrogen, neon) are too small to form 
either hydrates I or II because they ire not retained within the lat- 
tice very strongly. Byk and Fomina [10] have given a very extensive 
literature review on hydrates. 

Wilcox, Carson and Katz [11,12] proposed an empirical method for 
predicting hydrate formation conditions through the use of solid- 
vapor equilibrium ratios. This method supposes that the solid-vapor 
equilibrium ratio is analogous to the liquid-vapor equilibrium ratio 
and based on the idea that the hydrates form solid solutions. Further- 
more, the method assumes that the solid-vapor equilibrium ratio is 
independent of the overall composition of the mixture from which the 
hydrate has been formed. Even though it is now known that this 
premise is erroneous, the equilibrium ratio concept is useful in 
dealing with mixtures of reasonably uniform compositions such as those 
which may be found in natural gas production and transportation. Un- 
reliable results will be obtained whenever the hydrate-forming mixtures 
are greatly different from those for which the solid-vapour equilibrium 
ratios were obtained. A set of charts presenting equilibrium ratios 
for methane, ethane, propane, isobutane, carbon dioxide, and hydrogen 
sulfide is available in the Gas Processors Association Engineering Data 
Book [13]. The lack of experimental data on which the isobutane 
equilibrium ratios were based tends to make them subject to serious 
error. 

Van der Waals and Platteeuw [9] used statistical mechanical 


principles to derive the thermodynamic properties of gas hydrates 
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from a simple model by applying Lennard-Jones-Devonshire cell theory. 
McKoy and Sinanoglu [14] calculated and compared dissociation pressures 
of some gas hydrates using the Lennard-Jones 12-6, 28-7 and Kihara 
potentials in the Lennard-Jones-Devonshire cell model. They found 
that the Kihara potential predicted better dissociation pressures 
for hydrates formed from rod-shaped molecules. The LJ 28-7 potential 
model gives only a qualitative indication of the effect of narrowing 
the potential bowl on dissociation pressures, that is, in going from 
LJ 12-6 to LJ 28-7, the dissociation pressures become higher. The LJ 
12-6 potential does not depend on the shape and size of the interacting 
molecules. The model predicts good hydrate dissociation pressures 
for monatomic gases such as xenon and krypton but not for nonspherical 
and diatomic gases like CO, » No» 0, or CoH. 

Kobayashi and co-workers [15,16] modified the equations of van 
der Waals and Platteeuw [9] and proposed a method of calculating several 
thermodynamic properties of gas hydrates above the lower quadruple 
point by using LJ 12-6 potential. They obtained excellent agreement 
between the calculated and experimental data for pure gas hydrates. 
Later, Nagata and Kobayashi [17,18] used the Kihara potential to make 
Similar calculations for both pure gas hydrates and mixtures. For. 
this work they obtained parameters for the Kihara model from the litera- 
ture. Their results agree well with experimental data and are better 
than those obtained by using the LJ 12-6 potential. 

Parrish and Prausnitz [2] used a similar approach but they applied 
numerical methods to fit a new set of values of the thermodynamic 
properties of the empty hydrate lattice and the Kihara parameters for 


15 hydrate-forming gases. By using this method the agreement between 
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calculated and experimental data for many gas mixtures is better than 
it is for other previous methods. 

The hydrate formation condition for the binary system of methane 
and water has been reported in the literature by Villard CiSsz0}s 
Deaton and Frost [21,22], Roberts et al [23,24], Kobayashi and Katz 
[25,26], Campbell and McLeod [27], Otto and Robinson [28], and 
Marshall and co-workers [29]. Deaton and Frost [21] obtained the 
first reliable isobutane hydrate experimental data. Layano and Utida 
[30], and Rouher and Barduhn [31] have made a careful investigation 
of isobutane hydrate and determined its composition and quadruple 
points. Deaton and Frost [21] and Campbell and McLeod [27] investi- 
gated hydrate-forming conditions in the ternary system CHA-1-CQ Hy g- 
H,0 for a few low concentrations of isobutane. In all cases, the 


isobutane concentration was less than 4.6 mole percent. 
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CHAPTER 3 


THEORY 
A. Phase Behavior of Hydrates 

The Phase Rule is based on thermodynamic principles. It was 
first stated by Josiah Willard Gibbs [32] and it relates the number 
of phases (P), the number of components (C), and the number of degrees 
of freedom, or the variances (F) of a system in heterogeneous equili- 
brium by the mathematical relationship: 

FSG. = oP bo 

The only measured variables used in this equation are temperature, 
pressure and composition. The relationship deals with the heterogen- 
eous equilibrium encountered in processes which may be classified 
either as physical or chemical, as long as they are dynamic and re- 
versible in nature. In contrast it does not apply to the typical 
Static equilibrium states found in mechanics. - 

The Phase Rule deals with the phase behavior of a system as 
affected by changes in the variables for temperature, pressure, and 
composition. Composition refers to the proportions of the various 
components in the system and not really to the number and nature of 
these components. In a system where the composition is fixed over a 
certain region of pressure and temperature, the system behaves as if 
jt were a single component system. However, it does not follow that 
jt contains only one chemical species. 

In a two component system for example, consisting of one hydro- 
carbon and water, the phase rule states that with three-phase equili- 


brium such as hydrate-vapor-water rich liquid, there is only one 
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degree of freedom. If it is considered that one more hydrocarbon may 
be introduced into the above system in variable amounts for the same 
three-phase equilibrium then there would be two degrees of freedom. 
This would be represented as a surface on a phase diagram. On the 
other hand, if the hydrocarbon composition of this three component 
system is fixed, that is one phase yariable has been specified, then 
only the temperature or the pressure can be varied independently. 
Similarly, a multicomponent system in which the hydrocarbon composition 
is fixed has only one remaining degree of freedom for the three-phase 


equilibrium hydrate-water rich liquid-hydrocarbon rich liquid. 


Figure 1 is a qualitatively representative phase diagram for 
a system of three components with methane, isobutane and water. 
Methane is well above its critical conditions while isobutane is below 
its critical condition for three phases equilibria (H, Li» G). The 
line (Hy y> Ly> G)p is terminated by the quadruple point (Hyy> Ly» Los G), 
with the addition of a new phase, liquid hydrocarbon. Starting from 
the quadruple point, the four-phase equilibrium line is terminated by 
the critical locus of methane and isobutane because Ly becomes 
identical to G. This four-phase locus can be determined from the 
intersection of two three-phase equilibrium lines (Hyy> Ly; G) np and 


(Ly; L as indicated in the diagram or from direct experimental 


2° Sap 
measurements along the HiT Ly> Los G locus. 

Methane can exist in either hydrate Structure I or II, or both 
of them. When only hydrate I forming components are present in the 
system, methane will form hydrate with Structure I. But if there are 


some hydrate II forming components in the system, methane will form 
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FIG. 1 Qualitative Representation of Phase Behavior of the Ternary 


System CHy-1-CaHy g-Ho0 for Hydrate Formation 


* 


BTID -UIUHOAUO Stee 
Sel pT) io afd) 

ary dT eH) os 

bert pdt) fi 


tel slacevat 


: > ime 
au: et _,. ‘ 
i S 

a ‘Tr 


10 


hydrate II. Isobutane can only form in hydrate Structure II and 


consequently in mixtures of methane and isobutane only hydrate 


Structure II exists. 


B. Use of Statistical Mechanical Theory in Gas Hydrates Studies 


Van der Waals and Platteeuw [9] derived the solid solution theory 
of gas hydrates from a simple model which corresponds to the three 
dimensional generalization of ideal localized adsorption by using a 
grand partition function. They obtained several principal thermo- 
dynamic equations for gas hydrates such as the dissociation pressure 
of the hydrate at various temperatures, the heats of formation of 
hydrates, and the compositions of the co-existing equilibrium phases. 


The basic equations are: 
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where ‘od is the difference in chemical potentials of the empty hydrate 
lattice ie) and the filled hydrate lattice ui v5 is the number of 
cavities of hydrate type 7 per H,0 molecule in the lattice, Fa is the 
probability of finding a molecule of gas k in a cavity of a hydrate 
lattice of type 1; Ci is the Langmuir constant for molecules of k in 


cavities i, w(r) is spherical symmetrical potential of molecule k at 


the distance 


yosnt settufor rion oy ath hd: te] Ce sik efesh 7b. wey a 
bas ots 03 ebnogeay1o> fofdw fabon alamtz 6 mov? eatery e 

"Le antalya noi sorte bestinoof hasbi *o noissaiievensp’ heidi 
~uyrrants FegtanPaq favevae bsnfutdo yatT sac) toni oar Sem ie ~ 

a we29'Ng, atite tas ens 26 ive assebvd 26p VOT 2nof teuper atm 

0 mots 20 efor odd -zo1wdaysqmad enolic) J5 shorbyd ot - 


.eseaig tangy Betoine ss ont 70 2nottizognoa sit hn . come 
oS J 


i) 


ie 


a :976 cor yalps resi 


i 


CO Se Eee ae t4 OZ = tT) fats 5 TH- A ' Aj 7 Bia 
af A ae 


rt 


(7 a 7 +7) atu 3 TH 


i t 
See 46 eed be 6 es . Ase ) " = 
($) — tee ; + Pye 4) 
CB) ves netenn soni ABTA QW Ins, | = (i Yigg 


ote bya Wane sie Fo 2istinatog ‘writs ni sonsyotdhS aH ef hea : 
Yo "selvae ona et 7 Aa i) S31SIHf standby diel Ti SA bas (alas 
oid 2t gS « aabazet odd, nF sfuss tom Gall tag F days SBA! vo iP 
odarbyd & to wives 6 nt A 26p 46 Sfussvom 6 aria to Ny E> 
aa 7 oP, 
- oF 4 to hae seatrgh Ht ANNO. nett eit ect 63 Pa iter 
. ye . 
se fhm at cf aah ve 
UF ,) eee ; ; 


= 
of 
A 


Yr from its centre to the wall of the cavity, and ie is the fugacity of 
gas molecules j and is related to Y 59 the mole fraction in the gas phase, 


and p, the total pressure by 


Equation (1) is a generalization of Raoult's law which is for the 
properties of the solvent in a solution where solute-solvent interaction 
can be neglected. Equation (2) is Langmuir's isotherm for localized 


adsorption without interaction of the adsorbed molecules. 


Van der Waals and Platteeuw [9] considered the gas molecule moving 
about in a spherical cage formed by the water molecules and applied the 
Lennard-Jones-Devonshire cell theory. McKoy and S jnauog]u [14], Nagata 
and Kobayashi [17, 18], and then Parrish and Prausnitz [2] used Kihara 
potential for the force field in the cavity by consideration of the 
Shape and size of the interacting molecules. The L-J-D cell model 
with Kihara spherical core potential was modified as 
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where « is the potential minimum, o is the distance parameters where 


rr) 


r(r) = 0, and a is the core radius. In this potential model a 
molecule with a spherical core, e.g.; CH) » of radius a interacts with 
a “point molecule". The spherical core molecule represents the solute 
molecule in the cavity while the point molecule represents the water 


molecule in the lattice. 


Summing all of the gas-H,0 interactions in the cell, we obtain 
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where N = 4, 5, 10, or 11 in Equation 6, z is the coordination number 
and R is the cell radius of the cavity. | 

Using Equations (1), (2) and (6), we can calculate ne and then 
predict hydrate formation conditions, if we know the Kihara parameters 


€ and o. 


C. Numerical Method for the Prediction of Hydrate Formation Conditions 


For a hydrate to exist in equilibrium with the other coexisting 
phases the chemical potential of each phase should be equal. In the 
case ice is present, 


i debate (Tempe al! ©, J) PAVE dae tay pat. 7 (8) 


where u” (T, P) is the chemical potential of ice. If liquid water is 
present, 

H Sek 
fie hy oan ah rep eek en Xp otras esses eee eee e es (9) 
where we (T, P) is the chemical potential of pure liquid water at T and 


Pia ¢and Xv is mole fraction of water in the liquid phase. (x, is very 


close to unity for many gases.) When one lets Nr Oe ms = al. s0r 
au hs P) = ys ul and substitutes in Equation (1), then if ice is 
present 
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or if liquid water is present, 
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From Equation (10) or (11), we obtained the calculated potential difference. 
In this work fugacity coefficients were calculated by Chueh and Prausnitz's 
[34] computer program. Gas solubilities were estimated with the Krichevsky 
and Kasarnovsky equation [35] and the method of Morrison [36]. 

The chemical potential difference Au” or ie in Equations (10) or 
(11) respectively is a function of pressure and temperature only. When 
hydrate is in equilibrium with the other phases (ice and liquid water) 
the pressure and temperature effects on this chemical potential then 
become 

d(Au/RT) = = (AH/RTO)AT + (AV/RT)OP ...cccccescecveceseeee (12) 

where AH represents the difference between molar enthalpy and AV the 
difference between the molar volume of the empty hydrate lattice and 
liquid water. 

Integrating Equation (12) along the equilibrium line, it then 
becomes 


GREP URI Se Urea Ol RT = i" (AH/RT®)dT + E (AV/RT)+(dP/dT)dT ..(13) 


tis 1) 
where dP in Equation (12) is replaced by (dP/dT)dT because the system is 


univariant, and dP/dT is the slope of the pressure-temperature equilibrium 
curve. 

Using Equation (13), one can calculate the chemical potential 
difference [Au” (T, P) or Nis (T, P)] at the given temperature T and the 
reference hydrate dissociation pressure Pp for the reference hydrate when 
either ice or liquid water is present. 

If we let temperature be constant and integrate Equation (12) 
with respect to Pre we obtain, 
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The experimental chemical potential difference [au” (T, P) or 
ie (T, P)] can be calculated by using Equations (13) and (14). 

For temperature above 0°C methane and natural gas mixtures [22, 
27] were chosen as the reference components for Structure I and Structure 
II respectively. For temperature below 0°C, xenon and bromochlorodi- 
fluoromethane were chosen as the reference hydrates for Structure I and 
Structure II respectively. Pressure-temperature curves were fitted 
(2) for reference hydrates by the empirical equation 

Ine Pare oh 


+ Bp/T tr, alien 


R R R 
where Ap » Bp and Cp are constants fitted to represent the experimental 
data. These constants are given in Table 2. 

All the thermodynamic properties needed in Equation (13) and (14) 
are given in Table 3. AV was assumed to be independent of temperature 
and pressure. For the chemical potential difference between empty 
hydrate and ice at O°C and zero pressure, Parrish [2] used the reported 
value of 211 cal/mol [37] for Structure II. For Structure I, the 
previous value of 167 cal/mol [9] was doubted by Allen and Jeffrey [38] 
through their x-ray study of bromine hydrate because bromine hydrate has 
a tetragonal crystal structure rather than the cubic lattice of Structure 
I gas hydrates. A value of 302 cal/mol was estimated by Parrish [2]. 
This value compared well with the value reported by Child [39]. For 


the gases involved in this investigation, the Kihara parameters for 


gas-Ho0 interaction are listed in Table 4. 
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CHAPTER 4 


EXPERIMENTAL METHODS 


A. Experimental Apparatus and Materials 

(1) Equipment 

The equipment used in this investigation was essentially the 
same as that used by previous workers [28, 40] except for a newly 
designed trunnion used for rotating the cell and its contents. Figure 
2 shows a schematic diagram of the equipment and piping arrangement. 

The equilibrium cell A was designed and fabricated by the 
Department of Chemical Engineering and was similar to a Penberthy 
liquid level gauge. It was capable of withstanding 3,500 psia at 
100°F. The cell body was constructed from type 316 stainless steel 
and had an approximate volume of 80 cc. The front and back faces of 
the cell contained a pyrex window which was manufactured by the 
Corning Glass Company and was rated at 9,000 psia at 100°F. 

The cell A was enclosed in a lucite shel] B. The cell and shell 
were mounted on the trunnion C shown to the right of the main cell. 
A detail of the trunnion is given in Figure 3. The cooling fluid was 
circulated into the space between the lucite shell and the cell through 
the trunnion and back again to a temperature bath to control the cell 
temperature. The trunnion was connected to an agitation mechanism 
driven by a motor and acentric shaft which could be used to rotate the 
cell and its contents back and forth through an angle of approximately 
30 degrees. The cell was connected to the charging and sampling line 
on the top and to the tubing line going to the Ruska displacenent 


mercury pump P at the bottom of the cell. Mercury was used as the 
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confining fluid. The two tubing lines were made from spiral coils of 
Autoclave high pressure tubing inserted into the lucite shell at the 
left side. A SWAGELOK quick connection valve which was rated from 


vacuum to 2,500 psia was used at the sampling point. 


- 


(2) Control and Measurement of Temperature 


As indicated earlier, the temperature of the equilibrium cel] 
was controlled by circulating a coolant from an auxillary temperature 
bath through the lucite shell jacketing the cell. Denatured ethanol 
was used as the circulating fluid. The temperatures inside and outside 
the cell were measured by two copper-constantant thermocouples and a 
Leeds-Northrup Model 8686 millivolt potentiometer. The calibration of 


the thermocouples is shown below: 


Temp. (°C TUsErron: © 
Cel] Bath 
Ono =Caly77 -0.058 
24.4 -0.157 -0.109 
43.25 -0.115 -0.113 
97 .88 -0.403 ~0.445 


The temperature was controlled by an American Instrument Company bimetallic 
strip controller which could be set to govern either the refrigerator 
or the heater but not both simultaneously. By using this controller, 


the temperature inside the cell could be controlled to within +0.1°F. 


(3) Measurement of Pressure 
The pressure in the cell was generated by a Ruska Instrument 


Corporation positive displacement pump with a capacity of 250 cc and 
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rate at 10,000 psia. Two Heise Bourdon tube gauges rated at 

250 psia and 2,000 psia were used for pressure measurement. The 
gauges were calibrated with a dead weight tester. Corrections for the 
measured pressure in the cell were made for atmospheric pressure and 
for the head of mercury between the gauge and the cell. The accuracy 


of each gauge was within +0.05% of full scale. 


(4) Materials 

Matheson Incorporated ultra high grade methane and instrument 
grade propane and isobutane were used. They were certified to be 
99.97, 99.5, and 99.5 mole percent chemically pure respectively. The 
ethane used was Phillips research grade with a reported analysis of 
99.94 mole percent ethane. The analyses were confirmed by the fact 


that no detectable impurity peaks were observed during the chromato- 


graphic analyses of the samples. 


B. Experimental Techniques 
(1) Preparation of Mixtures 


Initially all the system was evacuated, flushed with the testing 
gas, and re-evacuated. Valve 7 was closed and approximately 10-15 cc 
of distilled water were introduced into the cell. Valves 11 and 6 
were closed, valve 10 was opened, and the predetermined amount of each 
gas was charged into the cell. The amount was based on the pressure 
reading observed on the proper Heise gauge. In making a mixture, the 
gas from the lower pressure tank was introduced first. With valve 5 


closed, the cell was continuously rocked for approximately 30 minutes 
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to ensure that the gases were thoroughly mixed. A sample of the mixture 
was then taken for analysis. 

When a condensed liquid sample was being prepared, it was found 
necessary to take several samples before consistent reproducible 
analyses were obtained. This was so because the contents of the lead 


lines could not be thoroughly mixed with the cell contents. 


(2) Hydrate Formation 


The cell temperature was lowered to 10-15°F below the estimated 
hydrated formation condition. More agitation of the cell caused the 
hydrate to form. When this had happened, the temperature was increased 
Slowly until the hydrate began to melt. The temperature was kept at 
0.3-3°F higher than the melting point to let almost all of the hydrate 
dissociate. The temperature was then lowered by 0.5-3°F to recrystallize 
the hydrate. A few crystals of hydrate that remained on the glass 
window served as seed crystals for the hydrate formation. The tempera- 
ture was raised again very slowly until the hydrate on the window just 
began to melt. Hydrate dissociation was always also accompanied by a 
sharp rise in pressure. The temperature and pressure were recorded at 
this point. The physical phenomenon of hydrate formation in the liquid 
phase resembles that in the gas phase. However, when the temperature 
was increased very slowly to the equilibrium point at a condition close 
to the bubble point in mixtures containing a relatively high concentra- 
tior) of methane, many bubbles developed in the liquid water phase no 


matter how carefully the procedure had been followed. 
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(3) Gas Chromatographic Analysis 


Sample analysis was accomplished by using a Hewlett-Packard 
Model 5750B gas chromatograph. Porapak Q columns were used with 
column length of three and six feet. The thermal conductivity cell 
was maintained at 120°C and 65°C, with 150 ma current and a Helium 
flow rate of 20 cc/min. The sample size was measured with a different- 
ial pressure transducer. The gas chromatograph was connected directly 
on-line with the Department IBM 1800 computer which permitted fast 
convenient data processing and digital filtering of data. The 
Chromatograph Monitoring Program was used to standardize each component 
investigated in the analysis to obtain the proper corresponding 
response factor, and to analyze the results of the samples. The 
smooth standardization curves are shown in Figure 4. By using the 
Chromatograph Monitoring Program analyses were reproducible to within 
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CHAPTER 5 


EXPERIMENTAL RESULTS 
A. Methane-Isobutane-Water System 

The experimental results for the different compositions of the 
methane-isobutane gas mixtures in the three-phase equilibrium (H, L, G) 
are given in Table 5. The eenults are plotted in Figure 5. In order 
to verify the experimental technique the binary system consisting of 
pure methane and water was tested. Good agreement between this work 
and that of other investigators [19, 20, 21, 22, 23, 24] was obtained. 
The experimental data of those authors and this investigation are given 
in Table 6 and depicted in Figure 5. For isobutane hydrate Rouher 
and Barduhn's data are tabulated in Table 7 and depicted in Figure 5. 
The experimental data of Deaton and Frost [21] and McLeod and Campbell 


[27] for several methane-isobutane mixtures are given in Table 8. 


B. hydrates of Liquefied Light Hydrocarbon 


Mixtures 

Initial hydrate formation conditions for the Ly-L,-H equilibrium 
for these liquefied hydrocarbon gas mixtures were determined. The 
composition of the mixtures and the experimental results are listed in 
Table 9 and Table 10. Smooth curves have been drawn through the 
experimental points as shown in Figure 6. The bubble point lines 
in the figure were predicted using a computer program based on the Chueh- 
Prausnitz correlation. It was an interesting and significant discovery 
to note that when the concentration of methane in the liquefied mixtures 


increased, the pressure had more effect on the hydrate formation. 
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temperature. Mixture I, containing no methane, had a nearly vertical 


equilibrium curve up to 1,000 psia. 


26 


mt : | i > a | > i pine 
, ! i at | i - ‘ os ibe 


Lia _ na : : 
. : Wty ) on; a 
7. a U9) Vesem ee 


rae 


CHAPTER 6 


THEORETICAL RESULTS 

A computer program based on the numerical procedure proposed by 
Parrish and Prausnitz [2] was written. The detailed procedure was 
described in Chapter 3, Section C. This program is given in Appendix C. 
At a later date, and through correspondence with Professor Prausnitz 
a copy of the program which had been used in the original paper was 
obtained. Both programs were used in carrying out the theoretical cal- 


culations. 


A. Hydrate Formation Prediction for Methane-Isobutane-Water System 


Predicted results compared with experimental data for methane and 
isobutane mixtures are given in Table 11. The program predicts well 
for higher mole percent isobutane mixtures but not so well for lower 
concentrations. No gas solubility data are available for isobutane in 
water in the literature. In preparing the program, it was assumed that 
the gas solubility of isobutane is approximately the same as that of 
n-butane. The gas solubility data obtained by Morrison [36] for methane 
and n-butane in water was used for the predictions of the methane-iso- 


butane mixtures. 


B. Hydrate Formation Prediction for Liquefied Hydrocarbon Gas Mixtures 
At the hydrocarbon bubble point, it should be possible to predict 


initial hydrate formation conditions in the liquid system because the 
vapour phase exists and is in equilibrium with liquid phase and hydrate 


phase. The hydrate formation equilibrium line in the liquefied phase 
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is extended to intersect the bubble point locus. The quadruple points 
(Ly; Los H,G) were located in Figure 6. The composition in the vapour 
phase at the bubble point was feed to predict hydrate formation condi- 
tions. A comparison between experimental and calculated hydrate forma- 
tion conditions is given in Table 12. The program predicted the hydrate 
fromation condition reasonably well with a maximum deviation in tempera- 
ture of -2.8°F. In a liquefied light hydrocarbon mixture containing a 
lower concentration of methane, the hydrate formation condition in the 


liquid mixture can be estimated without serious error. 
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CHAPTER 7 


DISCUSSIONS AND CONCLUSIONS 
A. Validty of Experimental Data 

Temperature measurements were accurate to within !0.1°F. Pressure 
could be determined to accuracies of +1.0 psi with the 2,000 psi Heise 
gauge and +0.125 psia with the 250 psi Heise gauge. By using the 
Chromatograph Monitoring Program, the sample analyses were determined 
with an accuracy with +0.3%. For the lower mole percent isobutane 
mixtures a Hewlett-Packard Model 17503A Recorder was used. The analyses 


were reproducible to within +0.02%. 


B. Methane-Isobutane-Water System 


From the theoretica] model we find the composition in the hydrate 
lattice is different from that in the gas phase. Therefore, the amount 
of hydrate remaining in the cell will effect the gas phase composition 
and consequently change the hydrate formation condition. It becomes 
very important to allow only a few hydrate crystals to exist when the 
point of incipient hydrate formation is recorded. In this way, the gas 
phase composition will keep approximately constant at different conditions 
for each mixture. For lower concentrations of isobutane, the effect of 
the amount of hydrate existing becomes even more significant. It was 
found necessary to analyse the vapour phase composition for each different 
condition. 

The data obtained by Deaton and Frost [21] deviate considerably 
from <his investigation. One of the reasons for the deviation might 


be that the equilibrium gas phase concentration of isobutane in Deaton and 
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Frost's work was actually somewhat less than that in the original 
mixture composition of 1.1 mole %. Another reason could possibly 


be errors in the analysis of the gas mixture. 


C. Liquid Hydrocarbon Mixtures 


From the experimental results it was found that methane has a 
Significant effect on the nature of hydrate formation for a liquid 
hydrocarbon mixture. In other words, in a mixture containing a relatively 
low concentration of methane, the hydrate equilibrium line in the liquid 
region will be nearly vertical. Therefore, if we can predict hydrate 
formation condition at the bubble point, we can estimate the hydrate 
formation at higher pressure for liquid mixtures without serious error. 
The experimental results confirmed the fact that the theoretical pre- 
diction method could also be used for liquid hydrate formation. A 
maximum temperature deviation of -2.8°F was found for a mixture contain- 


ing methane, ethane, propane and isobutane. 


D. Conclusion 

The study of hydrate of the methane-isobutane-water system has 
contributed to the available data on natural gas hydrate forming systems. 
Isobutane has been shown to exert a considerable influence on initial 
hydrate formation condition even when it is present in a concentration 
as low as 0.4 molez. 

The investigation of hydrate formation in liquefied light hydro- 
carbon mixtures verified that the theoretical prediction method could 
be used for representative mixtures containing variable amounts of 


methane, ethane, propane and iosbutane. 
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CHAPTER 8 


FUTURE WORK 


la 


The determination of the effect of carbon dioxide, hydrogen 
sulfide, and some of the non-hydrate forming heavier hydrocarbons 


on the L,-L,-H equilibrium. 


The completion of a computer program to predict hydrate 
formation at the bubble point for liquid mixtures by combining 
the hydrate formation prediction program with the bubble 


point prediction program. 
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NOMENCLATURE 
Core radius, A 


Fitted constants for calculating dissociation pressure 
of a reference hydrate 


Number of components 

Langmiur constant for components k in cavity 7 
Number of degrees of freedom 

Fugacity of component j 

Hydrate phase 

Hydrate phase with Structure I 

Hydrate phase with Structure II 

Molar enthalpy, cal/mole 

Gas phase 

Boltzmann's constant, 1.38 x 10°° erg/°K 
Liquid water rich phase 

Liquid hydrocarbon rich phase 

Number of hydrate forming molecules 


Number of hydrate forming molecules in smaller cavities 
of Structure I 


Number of hydrate forming molecules in larger cavities 
of Structure I 


Number of hydrate forming molecules in smaller cavities 
of Structure II 


Number of hydrate forming molecules in larger cavities 
of Structure II 


Integer constant in Equation (7) 
Number of "host" water molecules 


Total pressure, atm 
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Pe Dissociation pressure at ice point, atm 

Pp Dissociation pressure of reference hydrate, atm 
‘i Radial coordinate A 

R Gas constant, 1.987 cal/mole°K; cell radius, A 
T Temperature, °K 

ie Temperature at ice point, 273.15°K 

V Molar volume, cc/mole 

w(r) Spherically-symmetric cell potential, ergs 

a Mole fraction of liquid water 

ue Mole fraction of component j in gas phase 

Zn Coordination number of cavity type m 


Greek Letters 


T(r) Kihara potential, ergs 
5 Polynomial defined by Equation (7) 
é Depth of intermolecular potential well, erg 
ij Fraction of cavities type i occupied by component j 
m Chemical potential, cal/mole 
oF Number of cavities type i per water molecule in the hydrate 
O Distance parameter, point at which r(r) = 0, A 
ae Fugacity coefficient 
Superscripts 
) Property at ice point 


H Hydrate phase 
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Subscripts 


Liquid water phase 
Reference hydrate 
Ice phase 


Empty hydrate phase 


Component 


Cavity type i 
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APPENDIX A 
DATA PERTAINING TO HYDRATE PHYSICAL 
AND MOLECULAR STRUCTURE 
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TABLE 1 


PHYSCTAL PROPERTIES OF HYDRATE LATTICE (32) 


Structure I Structure II 
No. of water molecules/ 
unit cel] 46 135 
small Large Small Large 
No. of cavities/unit cell 2 6 16 8 
Cell diameter, A (9) 7.95 8.60 7.82 9.46 
Corrdination no. 20 24 20 28 
TABLE 2 


CONSTANTS FOR CALCULATING DISSOCIATION 
PRESSURES OF REFERENCE HYDRATES (2) 


(Pp is in atmospheres and T is in °K) 
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A -B C Teriperature 
_R oi! _R Range, °K _ 
Structure I 
Xenon 23.0439 -3357.57 -1.85000 211-273 
Methane -1212.2 44344 .0 187.719 273-300 
Structure II 
Bromochloro- 
difiuoromethane 11.5115 4092.37 O23 6033 253-273 
Natural Gas 
Mixture -1023.4 34984.3 159.923 273-291 


Natural Gas 


Mixture 4071.64 -193428.8 -599.755 eae Peco ls) 
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TABLE 3 


THERMODYNAMIC PROPERTIES OF EMPTY HYDRATE (p - Phase) 
- AND LIQUID WATER RELATIVE TO ICE (a - Phase) 


AT O°C AND ZERO PRESSURE (2) 


Structure I 


- 1 5 Cal/mole 
a Ni cal/mole 


= Vs cm?/mole 


= H”, cal/mole 


- Ch cal/mole 


TABLE 4 


302 


2f5 


onl) 


1436.3 


Ot 020336 iver eyo) 


Structure II 
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3.4 


KIHARA PARAMETERS FOR HYDRATE-GAS INTERACTIONS (2) 


Gas 


Methane 
Ethane 
Propane 
Isobutane 
Nitrogen 


co) 


Zoe n 


————— 


0.600 
0.800 
ie hen! 
1.600 
0.700 
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APPENDIX B 
EXPERIMENTAL AND PUBLISHED DATA 
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TABLE 5 


EXPERIMENTAL DATA FOR METHANE-ISOBUTANE MIXTURES - 
HYDRATE FORMATION IN HLG EQUILIBRIUM 


Composition (Mole % Of 


i-C,) psia 2p 

0.40 262 Shy 9) 
0.43 811 sok) 
0.46 457 46.0 
0.82 184 36.0 
ony 44.4 

485 5036 

855 57 26 

1,456 63.9 

Ihe i725 34.2 
201 40.2 

261 44. 

39] SOR 

503 Dome 

707 S/o 

1,008 (ne 

(ies) 101.6 34.3 
liye 40.3 

201 43.9 

one 50.3 

398 Sods 

50] Lye) 

662 60 

1,460 68.8 

6.0 VESTAS 34.9 
147.0 45.0 

245 52-5 

409 59.8 

Toe 44.) 3o 0 
81.8 42.3 

149.7 HOO 

295 60.3 

28.6 29.7 33.4 
SB ica 3952 

69.2 42.8 

87.3 45.7 

114.1 49.2 
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TABLE 5 - continued 
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TABLE 6 


COMPARISON OF METHANE HYDRATE FORMATION CONDITIONS (HLyG) 
OF THE PUBLISHED DATA AND THIS WORK 


Deaton and Roberts et al This 

Villard (19, 20 Frost,.0210u22 (23, 24) Work 
psia ad psia °F psia °F psia 
390 3220 4101S hss. 0 363° 32.0 408 
44) 34.0 421 34.0 848 46.0 509 
690 41.9 470 36.0 1,452 56.0 578 
934 47.3 496 37.0 1,567 56.4 668 
Tero? 49.8 553.) 39.0 713 
1,220 51.4 692 43.0 794 
1,813 ole 776 =©45.0 874 
5635 61.0 828 46.0 926 
2,650 63). | 879 47.0 15050 
3,410 66.7 982 49.0 1a 130 
3,900 68.5 halve. 3252.0 1,293 
1419 255.0 1,465 
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TABLE 7 


PUBLISHED EXPERIMENTAL DATA OF ISOBUTANE- 
HYDRATE FORMATION CONDITIONS 
IN HL,G EQUILIBRIUM (30) 


Sia we 
Toy SZ al 
17.9 32.8 
19.5 8000 
20.8 34.1] 
Cleo 34.7 


24.4 35.4 


44 


aie ter. willl 
iceblias Ld teas fr 


TABLE 8 


PUBLISHED EXPERIMENTAL DATA OF METHANE-ISOBUTANE 
HYDRATE FORMATION CONDITIONS 
IN HL,G EQUILIBRIUM 


Authors 


Deaton and 
Frost (20) 


McLeod and 
Campbell (27) 


Composition (mole %) 


Tal (Giniei -C 


1.4 
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TABLE 9 


COMPOSITION OF LIQUEFIED HYDROCARBON MIXTURES 


Nitrogen 
Methane 
Ethane 
Propane 


i-Butane 


mol %) 


TABLE 10 


EXPERIMENTAL HYDRATE FORMATION CONDITIONS FOR 
_____ FOR THREE LIQUEFIED MIXTURES 


Mixture I 
psia me 
Viz 40.2 
212 40.2 
324 40.2 
656 40.2 

1,026 40.3 


1,409 


40.4 


Mixture I] 
psig, ZF 
ia) 46.5 
357 46.6 
sie 46.9 
885 47.3 

1,194 47.8 
1,628 48.2 


Mixture III 
psia ae 
760 65.6 
856 66.0 
1,046 66.8 
1,421 68.0 
1,976 69.3 
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TABLE 11 


PREDICTED HYDRATE FORMATION CONDITIONS 
FOR METHANE-ISOBUTANE MIXTURES 
IN HL1G EQUILIBRIUM 


Composition Parrish and Prausnitz's 
mole % of Program (2) This Work Experiment 
1-C, ei psia psia psia 
0.40 Syaw) 323 oco 262.0 
0.43 ote heel 975 95] 811.0 
0. 46 46.0 544 542 UIE AO) 
O82 36.0 229 230 184.0 
44.4 400 399 eye) 
5OEG 616 609 485.0 
O70 sible! ige1g 855.0 
63.9 1 95 yale! 450.0 
hea oan 175.0 176.0 13765 
40.2 260 261 201 
Aen 338 Boe 261.0 
50.0 509 502 391.0 
yore 642 628 503.0 
57.8 923 883 707.0 
62-3 1,388 1,274 1,008.0 
pied Bas 129.1] 129.7 TQ 3:6 
40.3 19387 193.5 SV ae) 
43.9 247 246 201.0 
BOs 382 376 o1260 
537 487 A73 398.0 
SV eae: 636 608 501.0 
6180 858 797 662.0 
68.8 1,936 Wea 1,460.0 
6.0 34.9 90.4 90.5 MORE 
ASO 182.6 180.4 14720 
52.5 306 297 245.0 
59.8 519 483 409.0 
15.2 Sela(6) 52a Sysere 44.) 
A | 100.7 99.5 ea lets: 
50.5 181.9 176.0 149.7 
60.3 370 339 295.0 
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TABLE 11 - continued 


Composition Parrish and Prausnitz's 

mole % of Program (2) This Work Experiment 

1-Cy a psia psia psia 

28.6 33.4 SWE) o7,0 297 
39.2 60.0 59.4 on ley 
42.8 80.0 73.7 69.2 
noe 100.3 98.0 8753 
49.2 Ne lee Aster) 1421 

6226 33.2 23.0 23:./ 23.0 
See B22 31.9 $2.0 
38.8 41.6 41.0 Alpe 
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TABLE 12 


COMPARISON BETWEEN EXPERIMENTAL AND 
CALCULATED HYDRATE - FORMATION CONDITIONS 


AT BUBBLE POINTS FOR LIQUID MIXTURES 


Experimental Calculated Deviation of 
Liquid Mixture psia aa psia ig Temperature, °F 
I 169 40.4 157 3726 -2.8 
Ey 227 46.5 224 45.0 = hes) 


ita 667 bowl bog 63.0 =Z.4l 


- a : 
=p 4 : 7 77 

ns | a : . 
CMA HATHA AIANT WaawT as noaraiien 7 | 

oat Trias VOLT ara ~ SYAROYN GYTAdUI IAD 7 
2amTCIN ONIGE AOA TWIGS eeu A = 


75 eof fé ve a) 


APPENDIX C 
COMPUTER PROGRAM FOR GAS HYDRATE 
FORMATION PREDICTION 
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5] 
Flow Chart of Computer Program 


Read Pure Component 
Properties, 
T and all ue 


ae 


Estimate P 


Assume Hy or Hy 


Present 
t ie. 


Calc. Cy 
hin 


3{T 
P R) 


Set Pas P 


Calc . and x) 


P_.=P 
I] at T and P 


Assume Hy 


Present 
(10) or (11) for New P| 


| 
Solve Eq. (14) and oi | 


Hy 
Assumed 
Present 


| EES wit bast 
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MAIN PROGRAM 


THIS PROGRAM CAN BE USED TO PREDICT THE HYDRATE 
FORMATION CONDITIONS ABOVE ICE POINT FOR THE 
MULTICCMPOGNENT SYSTEM AND CALCULATE THE CGMPOSTION 
IN HYDRATE PHASE AT THAT CONDITION 


THE INPUT DATA STATEMENTS ARE DESCRIBED AS FOLLOWS 
NCOMP = NUMBER OF COMPONENTS 


KHC = CHOICE QF HYDRAE COMPOSITION COMPUTATION 
SET KHC=0 FOR GAS PHASE CALCULATION ONLY 
NTC = CHOICE OF FIXED TEMPERATURE 


SET NTC=0 FOR ONE TEMPERATURE CALCJLATIGON 


Di = INCREMENT OF TEMPERATURE WHEN NTC NeEew 0 
AL = NAME OF COMPONENT 

PC = CRITICAL PRESSURE 

vc = CRITICAL TEMPERATURE 

Ww = ACENTIC FACTOR 


CIRKV AND C2RKV = DIMENSTONLESS CONSTANTS IN REDLICH- 
KWONG EQUATION 

CORIJ = CHARACTERISTIC CONSTANT FOR I-J INTERACTIGN 

SAs SBs ANO SC = SGLUBILITY CONSTANT 

ALI, BLI, ALTI AND BLIii = LANGMUIR CONSTANT FOR 
HYDRATE STRUCTURE I AND STRUCTURE Ii 

Y = COMPCSITION OF GAS PHASE 


DIMENSION AL(455)s XM(4)sXS(4)5SA(4),SB(4)5SC(4) »S0(4) 
1 sHENRS(4)sHEN(4) »PMOLV(4) sALI (492) sALI1I( 492) sBLI(452)5 
2BL11(492)sCL1(452) »sCLI1I(452)sCL( 492) s TEMP (2) » TEMPY(4) » 
3SUM1 (2) sSUM2(2)UC2) sSY(4) sAN(4)5V(4) oF (4) sSCF( 2) sCF(4s 
42) +SML(492)HM(4) «UI (2) sUTI (C2) » PHF (4) 

COMMON PC(4)5VC(4) »TC(4S) oW(4) sCIRKV(4) sC2RKV(4) sTCIUV 
1(454) sPCIJUV(494) »CORIJ( 494) sNCOMPs PHI (4) 

READ(53100) NCOMPsKHCsNTCsDT 

WRITE(65100) NCOMPsKHCsNTCsDT 

DO 1 I=1i, NCOMP 

READ(55101) (ALC IsIN)»sIN=1;55) 

WRITE(63101) (ALC IsIN)»IN=135) 

READ(5,112) PCC(1I)sVC(1I)sTC(1) sWCT) sCIRKV(I),C2RKV(I)s 
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it 


333 
999 


1CGRIJ(141)sSA(1I)sSB(1).SC(1) 


WRITE (65114) PC(I),VC(1I)sTC(I) sW(1) sCiRKV(1I).,C2RKV(I), 


1CORIJ( 151) +sSA(1I)sSB(1I)3SC(1) 


IF (NCOMP-I1) 7.8537 

IIl=I+t+1 

DG 1 J=IIsNCOMP 

READ(5,104) CORIJ(I1I>J) 
WRITE(60104) CGRIJ( 1.44) 
CONTINUE 

CONTINUE 

DO 2 I=1,NCOMP 

Hoes t Jah. 2 

READ(53103) ALI(IsJ)sBLI(Is4) 
READ(52103) ALTII(1I,J),BLII(CI»sJ) 
WRITE(65103) ALICIs,/J)sBLI(CI.J)sALIICI2J)sBLII CI, 4) 
CONTINUE 

7O=273.15 

H=-2616.398 

DAH=20 .6166 

DBH=-C0 ..0211634 

DELVF=0.0387 

UI (1)=1./23. 

UEC2Si=Ss./e3< 

OLECTI=2e/ 47 

UII(2)=1./17. 

READ(5,102) TF 

IF (TF) 1121001511 

READ(S5,105) (YC 1) +I=1>+NCOMP) 
WRITE(62105) (YC I )sT=1sNCOMP) 
DG 333 I=1.NCGMP 
TEMPYCI)=YC(1) 

CONTINUE 

HI=0.0 

HII=0.0 

NOL=0 

ITER=1 

TR=TFt459.69 

TK=TR/1.8 


CALCULATE THE PARTIAL MOLAR VOLUME OF GASES IN THE WATER 


PH20=EXP(-2.4244+9.0363*(TR-492.)) 
PSOL2=14.696-PH20 

DO 31 I=1+NCOMP 

SY(1)=1.0 

DO 600 J=1,NCOMP 
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IF(J-1I) 610:;600,610 
eit SY CIl=0 50 
6CQ CONTINUE 
CALL PHIX(PSOL2Z2.TRsZVsSY) 
SOCI)=EXP(2.3026*(-SA(C I )#SB8(1)/TK+SC(1)*ALOG(TK)/S 
142.5026))) 
AN(CT)=SO(1)/(82.06*ZV¥TK) 
XMCT)=ANC(I)/55.51 
HENRS(1)=PSOL2*PHI(CI)/XM(I) 
PMOLV (I )=(0.095843. 76*TR*¥PC(1I)/(353156.*TC(1)))*(10.73 
1 9FGC1 Y7ECC Cr) ) 
31 CONTINUE 
DO 222 %t=1+sNCOMmP 
Y(I)=TEMPY (1) 
DO -222¢¢J=1 42 
CLICIeJ)=ALICIsJ)/TK#XEXPCBLI(I,J)/TK) 
CLIICIsJ)=ALIICIs J) /TKXEXP(BLIICI,J)/TK) 
222 CONTINUE 
c 
Cc TEST IF IT IS A PURE COMPONENT 
C 
DO 700 I=1sNCOMP 
MD=I 
IFCY(1)-1.0) 700,710,710 
700 CONTINUE 
GO TO 720 
710 IFC(CLI(MDs,2)) 73037403730 
720 HI==220 


GO TO 48 
740 HII=-2.0 
GGL TG. 201 
Cc 
2 TEST IF HYORATE I WILL FORM 
Cc 


720 DG 200 I=1,sNCOMP 
IF(CLI(I92)) 200,201,200 
200 CONTINUE 


HI=-2.0 
GO TO 48 
(e 
Cc 
C ASSUME HYDRATE I PRESENT IN WATER 
C 
48  PG=P 
HI=HI+1.0 


PII=P 
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DO Si I=1+sNCOMP 
“DO S51 J=152 
CL(IsJ)=CLI(I5J) 
CONTINUE 
RDELU=302. 
DELHI=275. 
DELHF=-1436.3 
DELVI=3.0 
DELVF=1.626 
AR=-1212.2 
BR=44344.0 
CR=187.719 
UC1)=UI(1) 
UC2)=UI (2) 

GO TO 30 


ASSUME HYDRATE [1 PRESENT IN WATER 


HII=HI1I+#1.0 

DO 4 I=1,NCGMP 
DON4 1 N= 52 

CL Cis JP=CLIICIs J) 
CONTINUE 
RDECU=S=2r1 . 
DELHI=193. 
DELHF=-1436.3 
DELVI=3.4 
DELVF=1.626 
UCHI=UFICH) 
U(2)=UTI(2) 

IF (TK—-291.) 59596 
AR=-1023.14 
BR=34984.3 
CR=159.923 

GO TO 30 
AR=4071.64 
BR=-193428.8 
CR=-599.755 


CALCULATE THE CHEMICAL POTENTIAL 
REFERENCE HYDRATE 


PGA=EXP (AR+ BR/TK+CR¥ALOG(TK )) 
PR=PGA 

CALL QG10(TGOsTKsARsBRsCRoFI ) 
FII=DELHI*(1.0/TK-1.0/TO) 
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40 


38 


39 


FII=FII+DAH*ALOG( TK/TO)+40BH*(TK—-TO) -H*¥(1.0/TK-1.0/TO) 
~PO=EXP( AR+BR/TOFCR¥ALOG(TQO)) 
RDELU=RDELU+( DELV I+DELVF) *P0*0.024214 
RDEL=TK*(RDELU/TO+F I14+(DELVItDELVF) *FI*0.024214) 
P=PGA*14.696 

PA=PGA 


CALCULATE GAS SOLUBILITY IN WATER 


EPS=0.01 

SUMXS=0.0 

NOL=NOL+41 

CALL PHIX(PsTRsZVsY) 

DO 32 I=1.NCOmpPp ) 
HEN( 1) =EXP (ALOGGC(HENRS(I1))+#PMOLV (1) *(P—-PH20)/(10.73* 
1 TR)) 

PHF (I )=PHI(1)*Y(1)*P 

F( I )=PHF(1)714.696 

XS(1I)=PHFE(I)/HEN(I) 

SUMXS=SUMXS+4#XS(1) 

CONTINUE 

XW=1.0-SUMXS 

DG 40 J=1.2 

TEMP(J)=0.0 

DO 40 1=1,NCOMmMP 
TEMP(J)=TEMPCI)FCLOCI5JS)*PHICI) *¥¥( I) 
CONTINUE 


USING NEWTON-RAPHSON METHOD TO SOLVE P 


ASUM1=0.0 

ASUM2=0.0 

SE 739 7J=132 

SUM1 (J) =U( J) X¥ALGG(1 -O+TEMP( J) *PA) 
ASUM1=ASUM1+4+SUMi (J) 

SUM2 (J) =U( J) X¥TEMP(C J) 01 .04+TEMP (J) ¥PA) 
ASUM2=ASUM2+SUM2( J) 

CONTINUE 
R1=1.987*TK*(ASUM1+ALOG( XW) )—-RDEL-(DELVI+DE_VF)*(PA-PR 
1 )*0.024214 
R2=1.987*TKXASUM2-(DELVI4tDELVF ) 0.024214 
R=R1/R2 

PA=PA—-R 

RA=R/PA 

ITER=ITER+1 

IFC ITER-200) 68026905690 
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RB=RA 

IF (LMT-200) 94361594 

IF CABS(RA)-EPS) 43,43,41 
IF (PA) 42542,38 

PA=0.5*PGA 

GO TO 38 

RB=(PA-PGA)/PA 

P=14,696*PA 

IF (ABS(RB)-0.001) 46546544 
PGA=PA 

IFC ITER-—200) 45,45,94 
EPS=ABS(0.1*RB) 

IF (EPS-0.01) 33,533,47 

EP S=O65.0 1 

GO TO 33 

IF(HII) 98:s49,49 

IF(CHI) 95350360 

LMT=200 

GO TO 48 

IF(P-PII) 97396561 

P=PII 

WRITE(65108) TFsPsITER»sNOL 
KIND=2 

GO TO 99 

RB=RBxP 

WRITE(63115) RB 

GO TO 45 

WRITE (60113) TFsPsITERsNGL 
KIND=1 

GO TO 99 

WRITE(63111) TF»Ps, ITER», NOL 
GO TO 99 

WRITE( 60113) TFsPsITERs NOL 
KIND=1 

GO TO 99 

WRITE(65108) TFsPsITER»sNOL 
KIND=2 

DG 170 I=1.NCOMP 

WRITE (69121) YCI) oC ALCIsIN) s IN=155) 
IF(KHC) 20935003209 

IF (NTC) 100031031000 
TF=TF-DT 

IF (TF-32.) 10,:9993999 
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209 IF(KIND-1) 205,205,207 

205 705206 °1=1.NCOMP 
DO 206 J=1,2 

206 ChOisJJ=CLI CIs) 
UC1)=UI C1) 
U(2)=UI (2) 
GO TO 240 

207 DO 208 I=1sNCOMP 
DO 208 J=1.2 

208° Ch OlsJ)=CLII Cis J) 
UC1jJ=UT1(1 ) 
UC2)=UI1G2) 
P=PII 

240; pO 210 7 5=1,2 
SCF(J)=1.0 
DO 210 I=1sNCOMP 
FCI )=PHIC1I) *Y(1)*P/14.696 
Ce Ciao) —cE( IoJ)*F( I) 
SCF (J)=SCF(J)F+CF( 154) 

210 CONTINUE 
SSM=0.0 
DO 212 1I=1.NCOMP 
SSML(1)=1.0E-8 
DG 211 J=1.2 
SML(IsJ)=CFCI,J)/SCF( J) 
SSML(1I)=SSMLC1I)4SML(CI5J5)*U(J) 

211 CONTINUE 

212 SSM=SSM+SSML(T) 
DO 220 I=1+NCGMP 

220 HM(I)=SSML(I)/SSM 
WRITE( 62120) 
DG 230 I=1+NCOMP 

230 WRITE(60121) HM(I)s(ALCI,IN)s IN=1,5) 
GO TO 500 

1061' GALE .EXET 

100 FORMAT(3155F10.0) 

101 FORMAT(SA4) 

102 FORMAT(F10.0) 

103 FORMAT(4E20.5) 

104 FORMAT(2F10.5) 

105 FORMAT(8F10.5) 


108 FORMAT(/5Xs*HYDRATE Ii WILt FORM AT T=*,F6.1" 
LF 8s 4Xs "P=" oF Bele" PSIA%37Xo ITH" s15s35Xo" OD=" 14) 
113 FORMAT(/5Xs*HYORATE IT WILL FORM AT =",F6.1i>5" 


1F* 54Xo "P=" sF8.15" PSIA® o37Xe*IT="* 9 1593Xs*'OD=" 514) 
sl piel FORMAT(/5X%9*BGOTH HYDRATE I AND II WILL FORM AT 
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LTS=*SFG6l4s 8 F654Xs *(P="SFB61s PSTA".20X%s"1T=".1553Xs 

2*OD=",15) 

112 FORMAT(7F10.0/3F10.0) 

114 FORMAT(7F10.5/3F10.5//) 

120 FORMAT(C//10Xe" THE HYDRATE CGMPCSITION IS*/) 

121 FORMAT (68XsF10.592X95A4) 

115 FORMAT(//10Xs" THE PROGRAM COULD NOT CONVERGE AFTER 200 
1 ITERATIONS® /15Xs*DEVIATON IS "sF9.25* PSIA*) 

END 


SUBROUTINE QG10 0 XL»sXUs ARs BReoCRsFI) 
FCT(T)=EXP (AR4+BR/T4CRXALGG(T))*(-BR/S T*X*¥34CR/TXX2) 
A= 6 5* (XU4XL) 

B=XU-XL 

C=.4869533xB 

Y= .03333567*(FCTCAtC) +FCT(A-C)) 
C=.4325317*B 
Y=Y4+.07472567* (FCT( AC) 47 CT (A-C)) 
C=.3397048*B 
Y=Y+.1095432*(FCT(A4C)4FCTCA-C)) 
C=.2166S977*8B 

Y=Y¥4+.1346334*(FCT(At+C) +FCT(A-C)) 
C=.07443717*B 
V=BX(V4214 77621 *(FCTCAtC) +FCT(A-C))) 
RETURN 

END 


SUBROUTINE PHIX(P5TsZVsY) 

DIMENSION ARKV(494) sBRKV(4) sAIRKV(4) sPHILN(4) sA(4) 5 
1 Z03) sY¥C4) sVCIIV06 454) sWI15 (454) ,ZC1I(454) 

COMMON (PC(4)sVC(4) sTC(4) »W(4) sCIRKV(4) sC2RKV(4) 5 
1TCIUV(494) sPCIUV( 454) sCORTII (494) »NCOMP sPHI( 4) 


CALCULATE VAPGR-PHASE FUGACITY COEFFICIENTS USING 
REVISED REDLICH AND KWONG EQUATION 


CALCULATE A AND B IN MODIFIED RK EQN FOR MIXTURE. 


DO 100 If=1,NCOMP 
ARKVCI,TI)=CIRKVCI )*10.73**2*( TCC) **2.5)/PC(T) 
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BRKV(CI)= C2RKV( I) *10.73*TC(I)/PC(1) 
IF (NCOMP -I) 111,5110.,111 
111 Il=frd 
DO 100 J=IIeNCOMP 
TCIIJVCTs JJ= (TCI) *TCOS) ) *¥*¥005*(1.0-CGRIU( 1.4) ) 
WISCIs J) =C WC IT) WS) )*O05 
2CTUCT sUJ=0.291=0-0S*wl JCI s J) 
VCIJV(Is J) =(VCC I) €VC(49) ) *0.5 
PEEIV Ot sy =2CIS Cis Ss ye1 0 2 7S tTCi IVC 563) SVCTIVCE. 5) 
ARKVCI sJ)=(CIRKVCI) tC1 RKV( J) ) *0.5%*10.73** 2xTCIUV(I oJ) 
1**2.5/PCIUV(IsJ) 
ARKV(JsI)=ARKV(I5J) 
100 CONTINUE 
110 CONTINUE 
AMRKV=0.0 
BMRKV=0.0 
DO 120 I=1,NCGMP 
AIRKV(1I)=0.0 
BMRKV=BMRKV4Y (1) *BRKV(I) 
DO 120 J=1,NCOMP 
AIRKV(I)=AIRKVC I) 4Y¥( JS) ¥ARKV(IsJ) 
120 AMRKV=AMRKV#Y(1)*Y( J) XARKV(IeJ5) 


CALCULATE VAPOR MOLAR VOLUME FOR MIXTURE, 


ACijJ=1.0 
A(2)=-1.0 
Ri 20.0 ed 
PBRT=P*BMRKV/SRT 

. ABRT=AMRKVZ (BMRKV*10.73*T**1.5) 
A(3)=PBRT*(ABRT-1.0-PBRT) 
A(4)=-ABRT*(PBRT*¥*2) 


CALL CUBEQ (AsMTYPE ;+Z) 


IF (MTYPE)130.1403,140 
130 IF (Z(1) -—Z(2)) 131.131,133 
131 IF €Z2(2)-— Z2(3)) 134.2134.132 
134° ZV=Z(C3}) 

60 TO°;15¢ 
132 zv=zZ(2) 

GO TO 150 
133 IF (Z¢€1)-Z(3)) 13421345135 
taS 2-764) 

GO TO 150 
Peon V=z2 CL) 
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150 VV=ZV*RT/P 
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CALCULATE FUGACITY COEFFICIENTS WITH MODIFIED RK 
EQUATION. 


QVVB=ALOG (VV/(VV=BMRKV) ) 

Q1VB=1.0/( VV-BMRKV ) 

Q@2RTB=2.0/(10.73*T**1 .5*BMRKV) 
QVBV=ALOG((VV+BMRKV)/VV) 
QARTB=AMRKV/ (10. 73% T*¥*1 oS *¥BMRKVXX2) 

QBVB=BMRKV/ (VV4tBMRKV) 

DO 160 I=1.,NCOMP 

PHILN(I )=QVVB+BRKV( 1) *QiVB-AIRKV( I) XQ2RTB*XOVBVt+BRKV(I) 


1 *QARTB* (QVBV-QBVB)-ALOG(ZV) 


1-QBVB)—-ALOG(ZV) 

PHI (1)=EXPCPHILN(1) ) 
RETURN 

END 


SUBROUTINE CUBEQ (AsMTYPE,Z) 


SOLVES CUBIC REDLICH=KWONG EQUATION FOR COMPRESSIBILI- 
TY FACTOR 


DIMENSION 8B(3)3sA(4).Z(3) 


BC1yY = (AC2Z)7AC1) 

B10V3 = B(1)73.0 

B(2) = AC3)7A(1) 

B(3) = A(4)7A(1) 

ALF = B(2) - 6(1)*B10V3 

BET = 2.0*B10V3**3 - B(2)*B10V3 + B(3) 


BETOV = BET/2.0 
ALFOV = ALF/3.0 
CUAOV = ALFOV *x*3 
SQBOV = BETOV **x2 


DEL = SQBOV + CUACV 
IF (DEL) 40,20,30 
MTYPE = 0 

GAM = SQRT (-ALFOV ) 
IF (BET) 22922321 

Z(1) = —-2.0*GAM -B10V3 
Z(2) = GAM -B10V3 

Z(3) = Z(2) 
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GG TO 50 

Z(1) = 2.0*GAM -B10V3 
Z(2) = -GAM -B10V3 
Z(3)0=-zZ@) 

GO TO 50 

MTYPE = 1 

EPS = SQRT (DEL) 

TAU = -BETOV 


RCU=TAUtFEPS 

SC U=TAU-EPS 

SIR=1.0 

STS=1.0 

IF (RCU) 31:32:32 

SIR=-1.0 

IF (SCU) 35434,34 

SIS=-1.0 
R=SIR*(SIR*RCU) **O0. 23333333 
S=SIS*(SIS*SCU) **0 . 32333333 
Z(1) = R + S = B10V3 


ZU2 a= —CRtS) 72.6 —-Bi0vV3S 
Z(3) = 0.86602540*(R-S) 
GO TO 50 

MIYPE = —1 


QUOT = SQBOV /SCUADV 
ROOT = SQRT (-QUOT) 
IF( BET) 42541541 


PEI = (1.5707963 + ATAN (ROOT 7 SQRT (1.0 - ROOT**2))) 
i 736S 

GO TO 43 

PEI = ATAN (SQRT (1.0 - ROOT**2) 7 RGOT) 7 3.9 


FACT = 2.0*SQRT (-ALFOV ) 


Z(1) = FACT*COS (PEI) - B10V3 
PEI=PEI42.0943951 

22) 2="PACTXCUS CPEL)} — B1TOVs 
PEI=PEI+4.1887902 

23) = FACT*COS CPEL) — S10VS 
CONTINUE 

RETURN 

END 
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